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ABSTRACT 



Context. Stellar evolution in close binary systems is strongly influenced by mass transfer from one star to the other when one com- 
ponent fills its zero- velocity surface or Roche Lobe. SSLep is a fairly nearby close binary showing the Algol paradox and a shell 
spectrum, both indicative of (past) mass transfer. 

Aims. To study the process of mass transfer and its evolutionary consequences, we aim at a direct characterisation of the spatial 
dimensions of the different components of SS Lep with IR interferometry. 

Methods. We use VINCI/VLTI interferometric observations in the K band and photometric observations from the UV to the far-IR. 
The visibilities are interpreted with simple geometrical models and the Spectral Energy Distribution (SED) is decomposed into the 
three main components: A star, M star and dust shell/disk. 

Results. From the SED, we find that the main emitters in the K band are the M star and the circumstellar environment. Both are 
spatially resolved with the VINCI observations, showing the excess to be circumbinary and showing the M star to have a size equal 
to its Roche Lobe. 

Conclusions. We conclude that we have, for the first time, directly resolved a star filling its Roche Lobe. The resulting mass transfer 
is probably the cause of (1) the circumbinary dust disk of which we see the hot inner region spatially resolved in our observations, (2) 
the unusually high luminosity of the A star and (3) the shell spectrum seen in the UV and optical spectra. 
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to binary motion is within the noise on the hipparcos measure- 
ments and that therefore the larger distance is the correct one. 



1. Introduction 

Close binary evolution differs from single star evolution in many 
respects. The most influential events after the main sequence 
(MS) are undoubtedly the mass-transfer episodes. These alter 
the evolution of both components, their chemical composition, 
their orbit and they often lead to spectacular high-energy phe- 
nomena. Needless to say, such mass transfer is a very complex 
event with many remaining questions regarding, amongst others, 
stability, mass conservation and flow shape. 

In her seminal paper. ICowlevI (Tl 967 ) suggests SS Leporis to 
be a semi-detached binary, consisting of a MS star and a giant 
which fills its Roche Lobe. The system shows the Algol para- 
dox, i.e. the most evolved component appears the least massive, 
which indicates past mass transfer. SS Lep is not strictly an Algol 
binary as it shows no eclipses. This is a consequence of its rather 
small orbital inclination. 

Before we discuss what is known about SS Lep to date, a 
preliminary note on the distance to this object is essential since 
the d ebate between Post- o r Pre-MS classifica t ion ha s been vivid 
(e.g. IWeltv & Wade1ll995h . |Poiidan & Shore! d 19911) mention an 
early trigonometric parallax measurement of 30 + 14 mas which 
places SS Lep as close as 33^ pc. However, the hipparcos par- 
allax is only 3.05 + 0.67 mas, corresponding to a distance of scenario of lCowlevId 19671) is therefore unlikely. 
330tgQ P c - Pourbaix (private comm.) confirms that any effect due 



Th e binary nature of SS Leporis was first detected bv lWilsonl 
(|1914|) in Lick Observatory spectr ograms. It w as found to have a 
variable absorption spectrum by Struve ( 1930) who also showed 
that the velocity variations of the sharp metallic lines are in- 
consistent with binary motion but rath er represent phenomena 
produced in an expanding atmosphere. Smith & Struve (1942) 
also found redward-displaced emission edges at the lo w Balmer 
memb ers and at the strongest Fen and Tin lines. ISlettebakl 
(Il950l) detected TiO lines in the IR spectru m of SS Lep in dicat- 
ing the presence of an M-type companion. Widing ( 1966) con- 
firmed the presence of this M compa nion and pro posed a Roche 
Lobe overflow model developed by iRrafj (Il958l) for T CrB as 
working hypothesis. 

ICowlevI dl967h was the first to determine the spectroscopic 
orbit (of the M star) accurately, see TableQ] Assuming an ad 
hoc radius for the M star of 75 R Q and an inclination of 24°, she 
found that the M star fills its Roche Lobe but only at per i astron . 
The orbital parameters were refined by IWeltv & Wadd ([l995) 
with a significant change in excentricity: 0.024+0.005 instead of 
0.132. The y concluded t hat the mass-transfer at/after periastron 
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Mid-IR emission features of silica te dust in t h e spec trum 
of SSLep were first identified by lAllen & Nevl d 19721) and 
lJura et all (1200 ll) suggested that SS Lep is surrounded by a cir- 
cumbinary dust disk with large grains, responsible for the far-IR 
and mm fluxes, from which an evaporation wind consisting of 
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Table 1. Relevant orbital parameters. P is the period, e the ex- 
centricity, K the velocity amplitude of the M star, a the semi- 
major axis of the Mstar orbit, i the inclination and f(M) = 
(M A sin0 3 /(M A + M M ) 2 . 





Cowlev (1967) 


Weltv& Wade (1995) 


P [days] 


260 ± 0.3 


260.34 ± 1.80 


e 


0.132 ±0.043 


0.024 ± 0.005 


K [km s 1 ] 


21.0 ±0.8 


21.32 ±0.21 


a sin i [10 6 km] 


74 


76.30 ± 0.53 


KM) 


0.24 


0.261 ± 0.005 



smaller particles arises and generates the excess at wavelengths 
below lOOyum. 

We report here on the direct size determination of the M star 
and the circumbinary matter using near-IR interferometry. 

2. Observations 

2.1. Photometry 

The photometric observations used to construct the SED were 
taken from the literature and span the entire wavelength range 
from the U V to the far-IR . Geneva 7-colour observations were 
taken from Rufener ( 1976 ). Near-I R magnitudes are those from 
the 2MASS PSC dSkrutskie et alj|2006l) and the far-IR magni - 
tudes were found in the IRAS PSC dNeugebauer et all 11984 . 
Additional observations in H,K and L are available fro m the 
CIT su rvey of early-type emission-line stars presented by lAllenl 
d!973l) . Variability in the optical is very s mall: 0.013 mag wit h 
a period of 130 days, i.e. exactly P or b/2 dKoen & Everll2002h . 
Variability in K is also very small with a detection at the 3 cr level, 
AK = 0.15 bv lKamath & Ashokl (fl999l 



2.2. Interferometry 

Interferometric observations in the near-IR (K band) were ob- 
tained with the VLTI commissioning instrument VINCI on base- 
lines with lengths ranging from 4.7 to 121 m and position an- 
gles (PAs) from 54 to 85 ° East-of-North. The 338 observations 
cover the period from 12 December 2002 to 30 December 2003. 
Instrumental visibilities were derived with t he standard data re- 
duction pipeline, v3.1 dKervella et al.l l2004). The visibility cal- 
ibration was done by linear interpolation of the interferomet- 
ric efficiency throughout the night. The stars used as calibrators 
are: v2CMa (K1III, 2.38 mas), 31 Ori (K5III, 3.56 mas), 51 Ori 
(Kim, 1.87 mas), HR2311 (K5III, 2.43 mas), HR2549 (K5III, 
2. 19 mas) and HR2305 (K3III, 1.76 mas). The Un iform Disk 
(UD) diameters were taken from Borde et al. (20021). 



3. Modelling 

In order to have some a priori knowledge on the angular sizes 
and the flux contributions of the different components in the K 
band, we first construct an SED. 



3.1. The SED 

We use the photometry presented in Sect. l2. II to construct the 
SED and fiQ a combination of a Kurucz model for the A star 



and a marcs model for the M star to the Geneva and 2MASS 
J,H observations. We searched for the best reproduction of the 
observations among the models (solar metallicity, spaced 250 K 
in r e ff and 0.5 dex in \ogg) having temperatures and surface 
gravities roughly compatible with t he spectral types (A l and 
between M3 and M4.5 respectively. IWeltv & Wadelll995l) . We 
chose to exclude the K an d L band observ ations in our fit 
since the SED presented by lJura et all (|2001|) and the VINCI 
observations (see Sect. 13. 2\ indicate some circumstellar excess 
already at these wavelengths. We do not fit the mid- to far- 
IR excess either since it is dominated b y cold dust emission 
(450 K, iFaiard o-Aco sta & Knacke 1995). The remaining free 
variables are the two stellar diameters and the total amount of 
extinction, i.e. both interstellar and circumstellar. The wave- 
le ngth dependence of t he extinction is modelled using the law 
of lCardelli etail(fT989l) with R v = 3.1. 

We find T eff = 9000 + 250 K and logg = 1.5 + 0.5 for 
the A star and T eff = 3250 ± 250 K and logg = 1.0 + 0.5 for 
the M star. We note that the Geneva colours agree much better 
with such a low surface gravity for the A star than with that of 
a MS star, which is independent evidence that the A star is in- 
deed located above the MS. A surface gravity of log #= 1 .48 was 
alread y derived from the Geneva colours by lHauck & Jaschekl 
(2000) but they neglected the influence of the companion. The 
M star's spectral type is rather uncertain, but we find a strong 
constraint on its temperature from our SED: the J and H band 
observations sample the peak of the intensity distribution (the 
H~ opacity minimum) in an M star and they are therefore very 
sensitive to its temperature. Our temperature is slightly below 
the temperature of a typ ical M4III gian t and corresponds more 
to spectral type M6III dPerrin et al.lll998l) . 

Our preferred model (x 2 =0.3) and its parameters are shown 
in Fig.Q] We find clear evidence for a near-IR excess in both the 
K and L band observations. In the Kband, at the wavelength of 
our interferometric observations, the A star contributes 11+ 7% 
of the flux, the Mstar roughly 60%. The remainder is unac- 
counted for in our model. 



3.2. The K band visibilities 

Several geometrical models are confronted with the data: a sin- 
gle UD, two and three concentric Uniform and/or Gaussian 
Disks and a binary model with possibly resolved components. 
A single UD yields a x 2 r — 39, and is thus clearly incompatible 
with the observations. A fit with two concentric UDs or a large 
Gaussian Disk and a smaller Uniform Disk, having three free pa- 
rameters: the 2 diameters and the flux ratio, has ax 2 = 6.4or5.5 
respectively. We should also include the unresolved contribu- 
tion of the A star. The baselines are not long enough to con- 
strain flux contribution and diameter of the A star so we use 
the parameters^ derived from the SED fit: #a.ud = 0.51 mas and 
F/t/F 'total — 0.1 1 ± 0.07 in K. This model, shown in Fig. [2] agrees 
rather well with the observations (x 2 = 5.5) and adding com- 
plexity, i.e. a possible offset between components, results in a 
degenerate problem: the x 2 hypersurface does not constrain the 
additional parameters (length and angle of the separation vec- 
tor), and shows an unreasonable minimal^ = 0.07. Fortunately, 
the derived diameters appear unsensitive to the introduction of 
such an offset. We must conclude that, while we can trust the 



1 The fitting is done by comparing observed integrated flux over the 
filter transmission profile with the same quantity computed from the 
synthetic SEDs. 



2 We computed the conversion factor from the Limb-Darkened (LD) 
diameter derived in SED fitting to the UD diameter used for the inter- 
ferometric modelling from our atmosphere models and found it to be 
0.945 in the K band for the M star parameters and 0.970 for the A star. 
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Fig. 1. SED, consisting of Geneva, 2MASS, CIT and IRAS pho- 
tometry, together with our model (solid line) consisting of the A 
and M stars (dotted and dashed lines respectively). Also shown 
is the blackbody curve at the brightness temperature (1250K) 
and diameter (10.5 mas FWHM) of the circumbinary dust disk 
generating the near-IR excess, as determined from the interfer- 
ometric observations (see Sect. 14.2b . To notice the excess in the 
Kband, one should realize that the bandpass actually covers the 
CO bandhead seen on the right of the datapoint. 

derived diameters, the current dataset does not allow for a char- 
acterisation of possible offsets (which are time-dependent since 
the orbital period is of the order of 1 yr), and this is therefore 
the aim of future AMBER and MIDI observations (see Sect.|5]). 
The inclusion of the barely resolved A star does not change the 
quality of the fit, but it does have a small influence on the de- 
rived diameter for the M star. We find the total uncertainty on 
the M star diameter and flux contribution with the inclusion of 
the uncertainty in A star flux to be 0.27 mas (10%) and 7% re- 
spectively. 

4. Discussion 

We summarize the results of our SED model and interferometric 
observations in Table|2] The flux ratios and diameters used in 
our geometrical model make an association with the components 
seen in the SED straightforward: the UD of 2.8 mas is the M star 
and the Gaussian Disk is the extended component generating the 
excess emission in K. 

4.1. The Mstar: filling its Roche Lobe 

Converting the measured UD diameter of the M star to a physical 
LD diameter yields # M LD = 3.11+0.32masorR M = 110+30R o . 
The uncertainty on the linear diameter includes that on the dis- 
tance. This is in good agreement with the LD diameter estimated 
from the SED fit. Given the absence of strong variability and 
the O-rich nature, we believe it to be on the Red Giant Branch 
(RGB). 

Th e constraint on the A star mass used by Weltv & Wade 
(I1995I) to derive the inclination should be relaxed somewhat 
towards the upper end since we have now the confirmation 



Fig. 2. VINCI/VLTI K band visibilities and our model (solid 
line) consisting of a Gaussian Disk and 2 Uniform Disks. The 
contribution to the total visibility of the extended component 
(Gaussian), the M star (UD) and the A star (UD) are shown with 
a dotted, dashed and dash-dot line respectively. Note that the 
A star diameter and flux and the uncertainties on these values 
were taken from the SED fit. 

Table 2. Stellar and CSE parameters derived from the SED 
and interferometric observations. NA indicates that this quan- 
tity could not be (reliably) measured, "idem" indicates that the 
SED-derived value was used. 



SED VINCI 



7rff,A [K], logg 


9000, 2.0 


NA 


r c ff,M [K], logg 


3250, 1.0 


NA 


a ,ld [mas] 


0.53 ± 0.02 


idem 


#m,ld [mas] 


2.66 ± 0.33 


3.11 ±0.32 


Shot excess [mas] 


NA 


10.5+0.1 


La [L q ] 


1900 ± 250 


NA 


Lm [Lq] 


1200 + 400 


NA 


F A at 2.2 /im 


1 1 ± 7% 


idem 


F M at 2.2 yum 


60 + 50% 


52 ± 7% 


F hot excess at 2.2 fim 


NA 


38 ± 7% 



that it is of luminosity class II. This yields an inclination of 
30° + 10°. Using (a A + a M )sin/ = 98.1 + 2.9million km or 
a = (141 + 4-Rq)/ sin/, and the generally assumed mass ratio 
Ma/ Mm = l/q = 4 ± 1, we find a Roche radiu^l ^R OC he,M = 
74^oR©, where the errorbar includes the 10° uncertainty on the 
inclination. 

From the observed surface gravity and radius of the M star, 
and the use of a reasonable mass for the A star (see Sect. l4.3t . we 
find that the mass ratio M A /M M could in fact be a little lower. 
This would reduce the size of the Roche Lobe for the M star. 

The find that the M star radius is equal to the Roche radius, 
which is very strong evidence that the M star is at this moment 
completely filling its critical Roche equipotential and therefore 
that mass transfer must be taking place. This is confirmed b y the 
shell spectrum with significant UV activity (Polidan 1994]), the 



3 **>°he = %JhCw*) ' << E Sgletoii|11983j) 
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presence of a circumbinary dust disk and the excess luminosity 
of the A star as shown below. 



4.2. The excess: a circumbinary dust disk 

The extended component we see in the interferometric observa- 
tions has a size clearly larger than the binary separation and we 
are therefore certain that almost half of the Kband flux comes 
from a circumbinary structure. We determine its brightness tem- 
perature to be Tb = 1250 K. This temperature is in agreement 
with the fact that we see the excess only from the near-IR on- 
ward and it corresponds more or less to the sublimation tem- 
perature of oxygen-based dust grains. This suggests that we are 
observing the hot inner region of the dusty circumbinary struc- 
ture also generating the mid- to far-IR excess seen in the SED. 
A clue to differentiate between spherical or disk-like geometry 
is the amount of reddening. We find Ay = 0.7 + 0.1 mag in our 
SED fit, w hich is in good agr eement with the E(B — V) = 0.26 
derived bv lMalfait et al.1 (1 19981) based on IUE data. This A v cor- 
responds to T,j ust ~ 0.5 at visual wavelengths. Since the line- 
of-sight extinction is either due to interstellar dust grains or due 
to circumstellar grains outside the orbital plane, we can assume 



th em to be small 
of iDorschner et al" 



We c an therefore use the optical constants 
(1 19951) for amorphous silicates to extrapolate 
the To.55^ m to IR wavelengths. We find an optical depth of only 
i"2.2/;m = 0.03 in the K band. The optical depth required to gener- 
ate the Kband excess at temperatures below that for dust subli- 
mation is of the order of unity, and we can thus conclude that the 
dust geometry is disk-like rather than spherical, and seen almost 
face-on, as expected from the system's orbital inclination. Note 
that a temperature much higher than the derived 7b would gen- 
erate excess emission also at shorter wavelengths, which is not 
seen in the SED. We can therefore rule-out the presence of a hot 
sp herical circumbina ry shell with significant continuum opacity. 

lJura et all (|2001) suggest a circumbinary dust disk with an 
evaporation wind to explain the mid en far-IR spectral properties 
of the excess. In their scenario the disk contains mainly large 
(d ~ 0.1 mm) grains. These would have an optically thin tem- 
perature at our derived distance from the central stars of 1480K 
which is in rough agreement with our derived value of 1250 K. 
The difference is probably due to the disk being fairly optically 
thick. The disk wind, consisting of smaller grains, is possibly 
responsible for the minimal amount of extinction still observed 
towards the central stars. 

The presence of this dusty disk shows that the mass transfer 
is not conservative which is in agreement with theoretical pre- 
dictions that the additional radiation pressure in the cool M star 
atmosphere may modify the critical Roche equipotential into a 
surface including both the inner and outer Lagrangian points 
dSchuermanll972l) . 

4.3. The A star: accreting 

The reaction of the primary, which has a radiative envelope, 
to the accr etion should be expansion, as its adiabatic exponent 
is positive (lRitterlll996l) . Indeed, we find a radius roughly ten 
times as large as that of a typical A1V star (Ra ~ 18 R Q ), 
while its mass, derived from the observed surface gravity and 
diameter, is not that of a (super)giant: .4M o < M A < 4M . 
Kippenhahn & Mever-Hofmeister d 1977b determined evolution- 



4 A spherical dust shell which is stable enough to allow significant 
grain growth is implausible: the radiation pressure would quickly drive 
them away. 



ary tracks for mass-accreting MS stars and find that the ex- 
cursion from the MS mass-luminosity relation depends on the 
accretion rate and initial mass. The position of the A star of 
SSLep in the HR diagram lies on the track of a 2M star ac- 
creting 2 x 10~ 4 M o /yr, after roughly 1 M has been accreted. 
For the mass ratio to become as large as is generally assumed 
{Ma/Mm = 4), this requires that the Mstar has also lost a sub- 
stantial amount of mass through the circumstellar disk or a wind: 
otherwise, the M star would have been the least massive to begin 
with, which is incompatible with its more evolved status. 



5. Conclusions and outlook 

We conclude that SS Lep is a semi-detached binary, consisting 
of a late-M giant which fills its Roche Lobe and therefore loses 
mass to the primary star. The latter was most probably a regu- 
lar MS star but now shows an increased size and luminosity due 
to a relatively high accretion rate. The mass transfer appears to 
be non-conservative since we detect an optically thick circumbi- 
nary dust disk. The shell spectrum detected at UV and optical 
wavelengths, in general not expected in A type stars, is thus most 
likely due to this ongoing mass transfer. We believe this to be the 
first time a Roche-Lobe filling star and the resulting circumbi- 
nary disk have been directly resolved using optical/IR interfer- 
ometry. 

We have planned further interferometric observations with 
AMBER and MIDI on the VLTI to characterize in more detail 
the shape and structure of both the M star and the circumstellar 
excess. We have also started a new spectroscopic campaign with 
the CORALIE instrument on the Swiss telescope at La Silla, 
Chile to improve the determination of the orbit and to study the 
gas flow dynamics. 

Algol-type binaries such as SS Lep can provide a wealth of 
information on mass transfer and binary evolution. Optical/IR 
interferometric observations such as those presented here, are 
well suited for studies of low-inclination systems and in that way 
complementary to eclipse and doppler-imaging work. 
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